Interplay of Magnetic and Superconducting Proximity Effects in FSF Trilayers 
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We present theoretical results on the interplay of magnetic and superconducting orders in diffu- 
sive ferromagnet-superconductor-ferromagnet trilayers. The induced triplet superconducting corre- 
lations throughout the trilayer lead to an induced spin magnetization. We include self-consistency 
of the order parameter in the superconducting layer at arbitrary temperatures, arbitrary interface 
transparency, and any relative orientation of the exchange fields in the two ferromagnets. We pro- 
pose to use the torque on the trilayer in an external magnetic field as a probe of the presence of 
triplet correlations in the superconducting phase. 



The importance of triplet pairing correlations in the 
interface region between a singlet superconductor and a 
ferromagnet recently became the focus of research in the 
field of spintronics [ll 12, S lj| ■ contrast to clean triplet 
p-wave superconductors and superfluids, for diffusive ma- 
terials p-wave correlations are suppressed and triplet cor- 
relations have s-wave orbital symmetry, but are odd in 
frequency . In the case of a homogeneous magnetiza- 
tion of the ferromagnet, the spin projection of the triplet 
correlations on the quantization axis of the exchange field 
is zero. If, on the other hand, the distribution of the 
exchange field in the ferromagnet is inhomogeneous in 
space, then under suitable conditions also triplet cor- 
relations with non-zero spin projection (equal spin pairs) 
are induced P, Q, 13 ■ Triplet pairing correlations induce 
in turn a spin magnetization both in the ferromagnet and 
in the superconductor • 

An important question is how to experimentally find 
good fingerprints of the triplet superconducting corre- 
lations. Theoretical work has been focused on calcula- 
tions of Tc H H 0, the local density of states (LDOS) 
0,0,[lllj_or to search for unconventional Josephson cou- 
plings |3, 1131 ■ Experimentally, no smoking gun has been 
found although recently magnetization changes were ob- 
served through neutron reflectometry on multilayers of 
YBa2Cu307 and La2/3Cai/3Mn03 (l3j. Various other 
properties of FS heterostructures that could be influenced 
by triplet superconducting correlations have been mea- 
sured, as for example a negative Josephson coupling (tt 
junctions) 0|, and LDOS modulations 

In the following we present results for the in- 
duced triplet correlations and corresponding changes 
of the magnetization in a ferromagnet-superconductor- 
ferromagnet (FSF) trilayer with arbitrary misalignment 
of the exchange fields in the two F layers. We put forward 
signatures of triplet correlations that can be measured 
experimentally. 

The FSF trilayer we consider is sketched in Fig. We 
denote the layer thicknesses by dpi, ds, and dF2 respec- 
tively. The X-axis is directed perpendicular to the layer 
interfaces with the origin at the center of the supercon- 
ductor. The z-axis is aligned with the exchange field Ji 
in the left ferromagnet. The angle between the exchange 



fields Ji and J2 is denoted Oj. We also assume transla- 
tional invariance in the y-z-plane. 

We use the quasiclassical theory of superconductivity 
for diffusive systems that is formulated in terms of 
momentum averaged Green functions. The Green func- 
tion for the trilayer, g{x, e„), depends on the x-coordinate 
and on the Matsubara frequency e„ = (2n -I- l)7rr (n in- 
teger, T temperature). In standard notation (see e.g. 
[3, llii ) the Green function is a 4 x 4 matrix in combined 
Nambu-Gor'kov (electron-hole) and spin space 



gs + 9f S 



ifs + ft ■ ^)i<Jy 



(1) 



where fs and ft are singlet and triplet pairing amplitudes, 
gs and gt are spin scalar and spin vector parts of the 
diagonal Green function, and the vector a = [ux, (Jy, (Jz) 
is composed of Pauli spin matrices. The Green function 
satisfy a Usadel-type equation 0, Hsf 



ienT3 



A- J- 



+ -dxigd^g) = 0, (2) 

TT 



where D denotes the diffusion constant (equal to Df in 
F and to Ds in S), [a, h] = ah — ba, and a = diag((T, a*). 
The quantities fi , 7^2 , tj, and 1 denote the Pauli matrices 
and the unit matrix in electron-hole space. Eq. ^ is 
supplemented by a normalization condition g^ = — tt^I. 
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FIG. 1: The trilayer consists of a superconductor of thickness 
ds and two ferromagnets of thicknesses dpi and dF2- The 
exchange fields of the ferromagnets, Ji and J2, are confined 
to the i/-z-plane, but are misaligned by an angle 9j. For 
torque measurements, a small magnetic field is applied in the 
y-z-plane at an angle ip. 
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The superconducting singlet order-parameter A is de- 
termined by the pairing interaction, which we assume to 
be zero in the ferromagnetic parts of the trilayer. By 
a proper gauge transformation, A can be chosen real, 
A — AicTj^fi, where A satisfies the gap equation 



(3) 



The pairing interaction and the frequency-sum cut-off 
have been eliminated as usual in favor of the bulk su- 
perconductor critical temperature Tco- The trilayer crit- 
ical temperature is lower < Tco, but we use Tqq as a 
parameter independent energy scale in our problem. 

The interface boundary conditions have been formu- 
lated by Nazarov 1^ and consist of two equations. The 
first is the condition of current conservation through the 
interface apgpdxgp = crsgsdxQs, where gp and gs de- 
note the values of the Green function at the F and S sides 
of the interface, respectively. The normal state conduc- 
tivities on either side are related to the corresponding 
diffusion constants D and densities of states at the Fermi 
level A// by cr = 2J\ffe^D, where e is the electron charge. 
Note that A// in the up and down spin bands are to qua- 
siclassical accuracy equal in the weak ferromagnet regime 
Tc < \J\ ^ Ef, where Ef is the Fermi energy. The sec- 
ond boundary condition is written as |2flj | 



o-pgpdxgF = ± 



1 



Stt^ [gp,gs] 



ARh 27r^2 ~T)^T{gF,gs}' 



(4) 



where Rb is the boundary resistance, A is the junction 
area, and T is the junction transparency. The sign +(— ) 
refers to the left (right) interface. For simplicity we use 
only two parameters: T and ri, = ARbcrp / (2Tr'^£^s), where 
^5 — ^ Ds /2ttTco is the coherence length in S. At the 
outer F surfaces the boundary conditions are dxgp — 0. 

We have solved Eqs. ©-lEl) self-consistently with the 
help of a parametrization of the Usadel Green function 
in terms of Riccati amplitudes, as described in Ref. [2ll |. 
This method allows us to treat spatially inhomogeneous 
(including non-coUinear) spin magnetizations. 

The magnetic order (exchange fields Ji, J2) and the 
superconducting order (A) are spatially separated. How- 
ever, due to quantum mechanical leakage of the supercon- 
ducting correlations trough the interfaces both the sin- 
glet fs and the triplet ft pair correlations are spread out 
through the whole structure. As a result also the magne- 
tization extends through the whole system, including the 
superconductor. The induced spin magnetization below 
Tr is defined as 



5M{x) = 2/^sAA/T^5t(x,e„), 



(5) 



where ^xb is the Bohr magneton. In Fig.[21we present the 
induced magnetization within the trilayer for parallel and 
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FIG. 2; Spatial dependence of the induced magnetization at 
various T G [0.05, O.SSjTco in steps of O.lTco for (a) parallel 
{9j — 0, Ji It J2) and (b) antiparallel {6j = tt, Ji ti J2) 
configurations. The S region is shaded, (c) Dependence of 
SMz on T at the left interface (in S) for both configurations, 
(d) Spatial dependence of the normalized LDOS in F2 at e = 
and at e = A(a; = 0) = l.l&Tco for the parallel configuration. 
Here, dpi = dF2 = ^s, ds = 4^s, Ji = J2 = 20Tco, ap = as, 
Dp — Ds, and both interfaces have rt = 0.1 and T — 1. 



antiparallel orientations of the exchange fields. We see in 
(a)-(b) that SM extends into the superconductor over 
the distance ^ £,s, while decaying and oscillating on the 
magnetic scale £,j — y/Dp/J in the ferromagnets. There 
is a large increase in the magnitude of SM as the temper- 
ature is lowered well below Tc, see Fig. ETc) . The corre- 
sponding spatial variations of the LDOS |22| are shown 
in Fig. |2Id). The LDOS is spin split as due to the pres- 
ence of triplet correlations. The order of magnitude of 
the LDOS modulations (1% effect) is in agreement with 
experiments |l5j . 

The singlet order parameter is suppressed at the in- 
terface to the ferromagnets. At the same time, the sin- 
glet correlations leak into the ferromagnets. For par- 
allel or antiparallel orientations of the exchange fields, 
triplet pairing correlations with zero spin projections are 
induced. For other orientations equal spin pairs are also 
induced. In Fig. |n{a)-(b) we show for perpendicular ori- 
entation the spatial dependence of the order parameter 
A{x) and of the correlation functions defined as |2^ 



n>0 

$tix) - /;(£„, x). 



(6) 
(7) 



n>0 



Note that because ft{—tn, x) — —ft{^n, x), the total sum 
over all frequencies in Eq. Q would vanish. Since we 
assume an energy-independent pairing interaction there 
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FIG. 3: Spatial dependence (a) of the singlet ($s) and triplet 
($t) pairing amplitudes, and (c) of the induced magnetization 
5M for perpendicular configuration (Ji _L J2, 6'j = 7r/2). 
Here, d^i = 8Cs, dF2 = dg = 4^5, and T = O.lTco. The 
other parameters are as in Fig. |21 The S region is shaded, 
(a) The singlet $s(a;) (dotted line) leaks into the F regions 
and oscillates and decays on a short length scale. The triplets 
^ty{x) (full line) and $tz(a;) (dashed line) are short-range or 
long-range depending on their respective projections on Ji 
or J2. Inset: spatial dependence of A(a;) in the S region, 
(c) There is an asymmetry in the oscillation period between 
5My (full line) and 5M^ (dashed line) in both F regions, (b) 
and (d): The angles 0$ and 9sm relative to the z-axis, that 
quantify the directions of and 5M in the ty-z-plane. 



is no triplet order parameter. The singlet component 
is purely real, while the triplet components are purely 
imaginary. The singlet component and the triplet com- 
ponent with zero spin projection on the local exchange 
field oscillate out of phase with respect to each other, and 
decay fast in the ferromagnets on the magnetic length 
scale "^^^ t)e seen in Fig. |3Ja)-(b), there 

are long-range triplet pairing correlations in both ferro- 
magnetic layers. The component ^ty is decaying slowly 
in Fi while ^tz is decaying slowly in F2. Relative to 
the local exchange fields, these components describe the 
equal spin pairing correlations that decay on the coher- 
ence length scale = \/Dp/2-kT but do not oscillate 
It should be noted, however, that these compo- 
nents are already at the interfaces quite small since they 
are induced non-locally in one ferromagnet by diffusion 
from the other through the superconductor. 

In Fig.|2Ic)-(d) we show the spatial dependence of the 
induced magnetization for perpendicular orientation of 
the exchange fields. Again, 5M is spread out over the 
large length scale ^5 in the superconductor. But 5M 
decays rapidly in the ferromagnet on the scale ^j, be- 
cause the spin-vector Green function gt is proportional 
to the fast decaying spin singlet component fs through 
the relation |^ gsgt = fsft- -^s a further consequence of 
this relation, the oscillation periods are different for the 



FIG. 4: (a) Change in magnetic moment 5m of a trilayer as 
function of the layer thickness dpi, for ds = 4^3, dF2 = O.l^s, 
T = O.lTc, and Ji = J2 = 20Tc. The other parameters 
are as in Fig. |5| (b) The temperature dependence of 5m for 
parameters as in (a), but with dpi ~ 0.2^s. The vertical lines 
in (a) are referred to in Fig. |^ 

two components of SM. For example, in Fi the oscilla- 
tion period of SMy is twice that of SM^. The reason is 
that 6My is determined by the product of the oscillat- 
ing fg and the monotonic fty in Fig.|2Ia), while (JM^ is 
determined by the product of the oscillating fs and the 
oscillating ftz- Note that deep in the ferromagnet (for 
X < — 5^s in Fi)j the magnetization change is mainly 
due to the long-range triplet correlations fty and SAI is 
therefore directed along the y-axis. But the magnitude 
of SAI is exponentially small in this region. An observa- 
tion of the two oscillation periods of the two components 
SMz and SMy, reflecting the different behavior of the 
short-range oscillating ft^ and the long-range monotonic 
fty, would be a smoking gun for long-range triplet com- 
ponents. Note that we have chosen a rather large value 
J = 20Tco in order to clearly separate the length scales in 
the problem. Smaller values of J might be more favorable 
in order to experimentally resolve this effect. 

Let us further address the issue of how to experimen- 
tally find fingerprints of the triplet superconducting cor- 
relations. We have seen that the local magnetization 
change is rather small. But it is important to realize 
that the integrated effect, 

Sm^Ajdx SM{x), (8) 

can be large. To illustrate this, we show in Fig. 0] the 
total magnetic moment of the trilayer, Sm, for various 
orientations of the exchange fields, both as a function of 
dpi and of T. We suggest to exploit the largeness of the 
integrated effect by measuring the torque on the trilayer 
in a weak external magnetic field. We consider a field 
B in the y-z-plane, aligned at an angle tp relative to the 
z-axis, as shown in Fig. ^b). Since the magnetic mo- 
ment m of the trilayer is also in the y-z-plane, the torque 
T = TO X B is directed along the x-axis. We assume that 
the external field is very weak so that the moment can be 
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FIG. 5: (a)-(c) The change in the torque relative to that in 
the normal state for parallel {6j — 0), antiparallel {9j = tt), 
and perpendicular {6j = 7r/2) orientations. The param- 
eters are the same as in Fig. 2] with the three numbered 
curves corresponding to the three layer thicknesses dpi — 
{0.01, 0.035, 0.21}^s (indicated by dashed vertical lines in 
Fig. |1J . In (d) we compare the total torque in the normal 
state (dashed line) with the one in the superconducting state 
(solid line from curve 2 in (c)) within the itinerant ferromag- 
net model. 



computed at i? = 0. In the normal state there is a torque 
due to the magnetic moment of the ferromagnetic layers. 
As the sample is cooled down through Tc, the magnetic 
moment changes, and consequently the torque is modi- 
fied. This is illustrated for the parallel and antiparallel 
exchange field orientations in Fig. |3Ia)-(b). In addition, 
when the exchange fields are not collinear, for example 
oriented perpendicular to each other as in Fig.|SJc), the 
(/?— dependence changes for asymmetric structures (here 
the layer thicknesses are different, dpi ^ c?f2)- The shift 
of the maxima in St^ {(f) when entering the superconduct- 
ing state, as e.g. in curve 2 in Fig. Oc), can be negative 
or positive depending on the parameters of the trilayer. 

In order to estimate the size of the change in the 
torque between the normal and superconducting state, 
we use the itinerant ferromagnet model. In this case 
the magnetization in the normal state is directly re- 
lated to the local exchange fields in the two layers as 
1/2- The corresponding magnetic mo- 



ments are — 2^B-^fJAdpi/2- The magnetic mo- 

ments in the normal and superconducting states are then 

5rh, respectively. 



'N 



'N 



and 



In Fig. IHd) we show how the torque is modified within 
this model. In particular, the equilibrium orientation (ip 
for which to || B) is different in the normal and super- 
conducting states and the maximum in Tx{ip) shifts to 
smaller (p by ~10%. The magnitude of these effects in- 
creases when parameters are such that the spin screening 
is more efficient (e.g. smaller J or smaller ^^1/2)- 



In summary, we have shown that superconducting 
triplet correlations in FSF trilayers are much easier de- 
tected by measurements of global properties of the sam- 
ple rather than by local probes. As an example we have 
studied the induced magnetic moment in the supercon- 
ducting state. The induced magnetic moment is a direct 
consequence of the presence of triplet pairing correlations 
in the trilayer. In an external magnetic field, the torque 
on the trilayer is modified in the superconducting state. 
We have shown that this effect can be more than 10% 
for exchange fields as large as |J| ^ 20Tco- The pro- 
posed effect can be used to experimentally determine the 
presence of triplet correlations. 
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